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Abstract. In this paper we present models for Single Stellar Populations (SSPs) of intermediate and old ages 
where dust enshrouded Asymptotic Giant Branch (AGB) stars are introduced. As long known AGB stars are 
surrounded by dust-rich shells of matter caused by their own stellar wind, which absorb the radiation coming 
^ ■ from the central object and re-emit it in the far infrared (IR). To this aim, particular care is devoted to follow 

the evolution of the AGB stars throughout the quiet and thermally pulsing regimes, to evaluate the effect of self 
contamination in the outermost layers by the third dredge-up mechanism, to follow the transition from oxygen-rich 
to carbon-rich objects (as appropriate to their initial mass and chemical composition), and finally to estimate the 
efficiency of mass-loss by stellar winds, all aspects that concur to the formation and properties of the dusty shells 
around. In addition to this, accurate physical models of the dusty shells are presented in which the re-processing 
of radiation from the central stars is calculated by solving the radiative transfer equations in presence of dust 
particles of different chemical composition. The resulting spectral energy distribution (SED) is examined to show 
how important features, like the 10/xm Si — O stretching mode feature and the ll/xm SiC feature, evolve with 
Oh 1 time. The SEDs are then convolved with the IRAS filters to obtain the flux in various pass-bands, i.e. 12, 25 and 

60 /im, for individual AGB stars of different, mass, chemical composition, and age. The comparison is made by 
5_l ■ means of SSPs along which AGB stars of the same age but different initial masses are located. This allows us 

to explore the whole range of masses and ages spanned by AGB stars. The theoretical results are compared to 
the observational data for selected groups of stars. The same is made for the J,H,K, L pass-bands of the Johnson 
system. Finally, from the integrated SEDs of the SSPs, we derive the integrated Johnson J,H,K, L magnitudes 
and colors to be compared to infrared data for star clusters of the Magellanic Clouds. In general good agreement 
with the data is possible if the effects of the circumstellar shells of dust are taken into account. 
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1. Introduction 

The first all-sky survey at far-IR wavelengths carried out 
in 1983 by the Infrared Astronomical Satellite (IRAS) 
has opened a new era of modern infrared astronomy. 
Thousands of galaxies were detected to emit most of their 
light in the IR. After IRAS, a long series of observations 
were started to explore the IR Universe. In the COBE 
all-sky maps a bright isotropic background (CIRB) was 
discovered in the far-IR/sub-mm spectral regions, whose 
origin is nowadays attributed to the integrated emission by 
dust in primeval galaxies absorbing and scatte ring the stel- 
lar lig ht and returning it at long wavelengths l|Puget et alJ 
Il996(l For the first time, the Infrared Space Observatory 
(ISO) surveyed and detected distant galaxies in the mid- 
and far-IR, allowing us to know in d etail their emission in 
those spectral regions ijElbaz et al.lll999|h 
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It became soon clear that precious information on the 
star formation history (SFH) of galaxies and the Universe 
as a whole is hidden in the UV-optical and IR ranges 
of the spectral energy distribution, thus spurring cross- 
correlated studies of the two regions. Many wide-field and 
all-sky surveys are currently running or have been just 
completed, e.g. the Galaxy Evolution Explorer [GALEX] 
llMartin et al.l ll997 ) and the Sloan Digital Sky Survey 
[SPSS] l|York et alJ 12000]) in the UV-op tical range; the 
Two Micron All-Sky Survey [2MASS] ()Skrutskie et alJ 
and the Deep Near-Infrared S urvey of the Southern 
Sky [DENIS] l|Epchtein et alJll997|i in the near infrared. 
Combined with other astronomical databases, they pro- 
vide a huge amount of UV-optical and near IR photo- 
metric data for millions of galaxies. In particular, the ob- 
servational information in the infrared will further grow 
with the new surveys SIRTF (Space Infrared Telescope 
Facility, see http://sirtf.caltech.edu/SSC/), UKIDSS (the 
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successor of 2MASS, see http://www.ukidss.org/) and 
with the advent of the Next Generation Space Telescope 
(NGST/JWST). 

This wealth of data must be accompanied by a con- 
tinuous upgrade of the basic theoretical tools to fully ex- 
ploit any information on the physical properties of the 
observed objects. In particular, to understand the role 
played by dust as strong IR emitter is mandatory. The 
key instruments for photometric studies are the SSPs, 
the building blocks of the assemblies of stars of different 
complexity going from star clusters to galaxies, and the 
fundamentals of population synthesis techniques. An SSP 
is defined as a coeval, chemically homogeneous assembly 
of stars all contributing to build up the integrated SED 
in a way proportional to the duration of their evolution- 
ary phase, luminosity and relative number (the luminosity 
function) , in other word s according to the precepts of Fuel 
Cons umption Theorem llTinslevlllQSfl lR enzini fc Buzzonl 
Il983l) . 

Let us shortly summarize here the various steps by 
which the theoretical SED of an SSP is derived, (a) given 
an age, the corresponding isochrone in the HRD is divided 
in elemental intervals small enough to assure that the lu- 
minosity, gravity, and T e // in them are nearly constant 
(the isochrone is approximated to series of virtual stars, 
for each of which we know the spectrum); (b) in each ele- 
mental interval the star mass spans a suitable range AM 
fixed by the evolutionary rate, therefore the number of 
stars per elemental interval is proportional to the integral 
of the initial mass function (IMF) over the range AM (the 
differential luminosity function); (c) finally, the contribu- 
tion to the flux at each wavelength of the spectrum by 
each elemental interval is weighed on the number of stars 
in it and their luminosity. To conclude, basic ingredients 
of the SSP SEDs are the isochrones and their path in HRD 
(in turn functions of the initial chemical composition) , the 
IMF, and a library of stellar spectra for different values of 
T e ff, gravity, and chemical composition. 

With a few exceptions to be mentioned later, the light 
emitted by SSPs is m odeled neglecting the presence of 
dust around their s tars (jBertelli et alll994llTantakTl 998: 
iGirardi et al.ll2f)f)2l e.g.). However, there are at least two 
situations in which dust is present: the very initial stages 
when stars are born embedded in molecular clouds, and 
the late stages of AGB stars when dusty envelopes are 
formed. In this study we will concentrate on the AGB 
stars leaving to a forthcoming paper the case of young 
stars. 

The dust shells surrounding AGB stars are the result 
of mass-loss by stellar winds and the complex structure 
and evolution of these stars during the thermally pulsing 
AGB phase (TP- AGB). They trap the radiation coming 
fro m the central AGB star and re-emit it in the far IR 
rsee lHabinglll99fil for a classical review of this topic). The 
contribution to the IR radiation by the dusty shells of 
AGB stars is particularly relevant because AGB stars are 
luminous objects able to significantly affect the integrated 
SED of star clusters and galaxies. 



A great deal of studies have been devoted to model 
the circumstellar shells around AGB stars and to under- 
stand in detail important properties of the shell structure, 
such as the composition of the dust grains, the dust-to-gas 
ratio, the exp ansion velocit y of the matter, and the mass- 
loss rates fcf. lHabin"gl [l996'). However, only an handful of 
studies have tried to include the effect of the dust shells 
surrounding AGB star s on the integrated S EDs of SSPs: 
the pioneer paper by IBressan et al.1 lll998h and the re - 
cent series of three ar t icles b y IMou hcine fc LanconI l|2f)02h : 
lLancon fc Mouhcinel l|2002|) and iMouhcinel (|2002|) . 

There are several causes hampering detailed studies of 
the subject: 

(i) First of all, theoretical spectra of oxygen-rich (O- 
stars) and carbon-rich (C-stars) AGB stars, surrounded 
by shells of matter and even of simple AGB stars are not 
available. Only empirical spectral are to our disposal. This 
is an obvious difficulty for any fully theoretical analysis of 
the problem. The spectra of AGB stars of any type are 
particularly difficult to obtain because of the many pa- 
rameters entering the problem, i.e. the spectrum of the 
inner stellar source, the optical depth (key parameter of 
the radiative transfer problem to be solved to follow the 
propagation of radiation throughout the dusty shell) , and 
the properties of dust grains in the shell. These latter de- 
pend on the optical depth in such a way that we have 
different compositions for different optical depths and of 
course passing from O-stars to C-stars. 

(ii) Second, for long time SSPs have been developed 
to interpret observational data in the optical range of the 
spectrum. Only recently the wealthy of data in the IR 
have spurred several groups to develop sui table SSPs ex- 
tendin g to the far infrared. The study of IBressan et alJ 
(I998J, despite severe limitations such as the use of gi- 
ant spectra to model the radiative transport in O-stars 
and C-stars, the lack of the transition between O-stars 
and C-stars, and the crude treatment of the dust compo- 
sition, has op ened the way to a new ge neration of SSPs 
in the far IR. IMouhcine fc LanconI l)2002f) and IMouhcinel 
(2002) have improved the situation including appropri- 
ate models for the AGB phase, the transition between 
O-stars and C-stars, and an empirical library of spec- 
tra of long period variables (LPV) as input for the cir- 
cums tellar shells. I t is im portant to remind here that 
while IBressan et all lll998l) made us e of purely theoretica l 
spectra. IMouhcine fc LanconI (|2002h and lMouhcinel l|2002h 
have adopted an empirical library of spectra of AGB stars, 
so that part of their results is implicit in this initial as- 
sumption. 

Aim of this study is to go further along the theoretical 
line of work and to generate modern integrated SEDs of 
SSPs that: (i) extend to the far IR; (ii) include more up- 
dated models of AGB stars; (iii) allow for the metallicity 
dependence and finally (iv) include an accurate modelling 
of the shells of dust surrounding these stars, thus improv- 
ing upon the emission in the IR spectral r egion. We will 
follow a purely theoretical approach as in IBressan et alJ 
(1998) and will try to include a detailed treatment of the 
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circumstellar dusty shells and of the transition from O- 
to C-stars. It is clear that the theoretical spectra of M gi- 
ant stars are not as good as the empirical ones (they can 
indeed be applied only to the case of unobscured O- and 
C-stars). However, we will insist on the purely theoretical 
approach in order to establish how far current theory can 
go in interpreting and reproducing the observational data. 



This paper is organized in the following way. In Sect. [21 
we present the model for the envelope of AGB stars, de- 
scribe the assumptions for the radiative transfer problem, 
and calculate an expression for the optical depth, the key 
parameter of the radiative transfer equations, and its de- 
pendence on basic stellar parameters such as the star mass 
M, radius R, effective temperature T e ff, luminosity L, 
pulsation period P, and metal content Z. In Sect. [3| we 
briefly summarize the prescription we have adopted for 
the mass-loss rates of AGB stars. In Sect. 0J we discuss 
in detail the optical properties of dusty AGB shells. As a 
matter of fact, the dust absorption coefficients are differ- 
ent for O-stars and C-stars, and even for the same type of 
stars, different optical depths may imply different proper- 
ties. In Sect. we present our choice for the parameters 
governing the dusty shells of AGB stars, i.e, the temper- 
ature on the inner boundary, the density profile across 
the shells, and the type and mixture of grains. In Sect.[Hl 
we try to discuss and evaluate the uncertainty affecting 
the models of the dust shells and how this would reflect 
onto the SED of the SSPs with particular attention to 
the IR range. In Sect.[7|we define the monochromatic flux 
emitted by a SSP, summarize the assumption for the ini- 
tial mass function and the rates of mass-loss from stars 
that we have adopted, report on the libraries of stellar 
models, isochrones and stellar spectra with aid of which 
we have calculated our SSPs, and finally discuss in detail 
the problem of the transition between O-stars and C-stars 
AGB stars and how this is incorporated in our SSPs. In 
particular we describe how we have included the recent 
models b y iMarigo et al.1 (jl99flft into the SSP models by 
iTantalol l|l998h . In Sect. we present the SEDs of our 
new SSPs that take the effect of dusty AGB shells into 
account an d compare our results with those for the stan- 
dard SSPs l(Tantalolll99^ . In Sect. El we describe the in- 
tegrated IRAS far IR colors predicted by our models and 
compare them with the observational data for a sample 
of AGB stars (Mira, Semi-Regular Variables, Long Period 
Variables, OH/IR stars and C-stars). We also derive the 
integrated near IR colors of the SSPs and compare them 
with the data for a selected sample of young star clusters of 
the Magellanic Clouds. In addition to this, we examine the 
age dependence of the inte grated near IR color s of SSPs 
from different sources, i.e. iBerte lli et al .1 lll994l); Tantalol 



Jl998l):lGirardi et all feOOak iMouhcine fc Lanconl l)2002|) ; 
lLancon fc Mouhcinel l|2002h and this study. Finally, some 
concluding remarks are drawn in Sect. ITU1 



2. Modelling a dusty envelope 

The problem of radiative transfer in the dusty shells sur- 
rounding the mass losing stars (AGB stars in particular) 
has been addressed by many authors (see lHabindri996l 
and references therein). The most complete formulation 
of the problem couples the radiative transfer and hydro- 
dynamical equations for the motion of the two interacting 
fluids, gas and dust, and takes the interdependence be- 
tween gas, dust and radiation pressure into account. 

Since our goal is to simply include the effects of the 
AGB star dusty shells on the spectra of SSPs, we will limit 
ourselves to consider the problem of radiative transfer in 
the shell and will leave hydrodynamics aside. To this aim, 
we need to link the fundamental parameters of a star, i.e. 
mass M, radius R, luminosity L, pulsational period P, and 
metal abundance Z to the parameters that characterize 
the shell of matter and that are relevant to the solution of 
the radiative transfer across it. 

The radiative transfer equation is 



dh (I) 
dl 

where 



= k x (0 [s x (0 - h (i)} 



(1) 



k\ (I) is the overall extinction coefficient at the wave- 
length A given by the sum of the absorption and scat- 
tering coefficients, kx (I) = k a x (I) + k s x (I); 
Ix (I) is the intensity of the radiation field; 
Sx (I) is the "source function" , given by the ratio 
£\ (0 /kx (I), with ex (0 the emission coefficient. 



As demonst r ated bv | Rowan- Robinson! (^980) and 
llvezic fc Elitzurl (|l995l ri997[) 7 only two dimensional scales 
are involved in the radiative transfer equation because Ix 
and Sx have the dimension of an intensity, whereas kx 
and I have the dimension of a length (kx is the inverse of 
a length). Therefore, any physical quantity related to the 
radiative transfer problem can be expressed as functions 
of these two scales. 

Defining the dimensionless element of optical depth 
along the ray path dl as drx = kx (/) dl, we can write eq. 
(Q) as 



dh (I) 
dr x 



Sx (0 - h (0 



(2) 



Starting f ro m these simple considerations, 
llvezic fc Elitzurl <|l997l) have proved that the radia- 
tive transfer equation satisfies the property of scale 
invariance: the physical dimension of any system can 
be increased and decreased in an arbitrary way without 
affecting the radiative properties, as long as optical 
depths and spatial variation of the opacity remain 
the same. Two systems with different dimensions and 
absorption coefficients, but with the same total optical 
depths and auto-similar distributions of opacities and 
"source functions" will produce the same intensity of the 
radiation field. 
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iRowan-Robinsor] ljl980|) first applied the radiative 
transfer scale invariance to the 1R emission of a cen- 
tral source surrounded by a dust shell. Subsequently, 
llvezic fc Elitzurl l|l997h presented a general formulation 
of the problem in arbitrary geometry and distribution of 
dust, and studied in detail the case of a spherical shell of 
dust heated up by a central source. They also pointed out 
that the concept of scale invariance is particularly use- 
ful when the absorption coefficient does not depend on 
the radiation intensity. Unfortunately, as a c onsequence 
of this, the analysis bv llvezic fc Elitzurl cannot be 

applied to emission or photoionization lines (where the 
absorption coefficient can depend on the intensity of the 
radiation field via its effects on the level populations) but 
only to the continuum of the radiation coming from dust 
heated by a central source. Last point to note is that when 
the shell is optically thin over all wavelengths A, the ap- 
proximation t\ ~ may be used thus obtaining a quick 
analytical solution of the problem. In contrast, when the 
shell is optically thick the problem must be solved numer- 
ically. 

In the following we apply t he method, the results a nd 
the numerical code DUSTY bv llvezic fc Elitzurl l)l997|) to 
study the dust shells surrounding AGB stars. The spher- 
ical symmetry approximation is adopted for the sake of 
simplicity. The key parameter we need to solve the radia- 
tive transfer problem and to calculate the emerging flux is 
the optical depth T\ of the shell which is defined as follows 



dr x (r) = I k\ (r) p (r) dr 



(3) 



where k\ is the overall extinction coefficient per mass unit, 
p is the matter density. Both depend on the radial distance 
r from the central source. The integral is evaluated over 
the thickness S of the shell. In the case of a dusty shell 
eq. (PJ) becomes 



T x = k x (r) p d (r) dr 
Js 

where pd (r) is the dust density. 

The equation of mass conservation is given by 



(4) 



dM (r) 
dT~ 



47tr 2 p (r) v (r) 



(5) 



The matter density p (r) is linked to the dust density pd (r) 
by the relation pd (r) = p(r) S where S is the dust-to- 
gas ratio. Substituting the matter density p with the dust 
density pd we obtain 



dM (r) _ Airr 2 p d (r) v (r) 
dt ~ 5 

and then 



Pd 0) 



M (r) 5 
4irr 2 v (r) 



(6) 



(7) 



With aid of the above relation we recast the optical depth 
of eq. as 



r A = / fc A (r) 
Js 



M (r) 5 
4irr 2 v (r) 



dr 



(8) 



In reality, the dusty shell will extend between an 
internal and external radius, rj n and r out respectively. 
Assuming 8 to be constant across the shell, we get 



or 6 r° Ut , , , MM , 

T A - — / k x (r) ^y,dr 
Air J r r z v (r) 



(9) 



To proceed further, the mass-loss rate M(r), velocity 
v(r), and extinction coefficient k\(r) and their radial de- 
pendence must be specified. To a first approximation we 
assume that at any given time the rate of mass-loss and the 
velocity are constant with r, together with physical prop- 
erties of the dust. Therefore the absorption coefficient for 
unit mass kx does not depend on the radial coordinate r. 
With these simplifications we have 



T x = 



sMkx r° ut i 



-rrdr 



Attv 

and upon integration 

_ SMkx ( 1 I 
1\ 



Airv V r ln r out 



(10) 



(11) 



Since r out is usually much larger than ri n , the relation 
above can be safely approximated to 



SMkx 1 
Attv r in 



(12) 



Now we need to connect the physical quantities defin- 
ing the optical depth of the dusty shell with the typical 
parameters of AGB stars, e.g. mass M, radius R, effec- 
tive temperature T e f f , period P of pulsation, luminosity 
L, and metal content Z. Assuming the shell to be opti- 
cally thick to IR radiation (a good approximation in case 
of high mass-loss rates), the inner radius of the shell can 
be derived from the equality 



ft 



47T7-L CrT? 



(13) 



where i?* is the stellar radius and Td is the temperature 
of dust condensation at fj„. From this relation we get 



( L, 



L 



(14) 



The uncertainty arising from using the relation l|13fl 
can transferred to Td which is not firmly estab- 
lished. Literature values r ange from 800 to 1500 



1! 



Rowan-Robinson fc Harris! Il982t iDavid fc Pap oularl 



1990; Suh 1999, 2000; lLorenz-Martins fc Pompei 
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lLorenz-Martins et al-lkoOll: ISuhll2002|) . Finally, adopting 
T d = 1000 K (see Sect. below) we obtain 



= 2.37-10 



12 



L 



(15) 



The extinction coefficient for unit mass k\ is in general 
given by 



Pd 

where the summation is extended over all types of grain 
in the mixture, a gi is the cross section for the i-th type 
of dust grain and n 9i is the number of grains of i-th type 
for unit volume. Denoting with m Bi the mass of the i-th 
type of dust grains, and introducing the mass abundance 
Xi — n gi m gil Pd of i-th type of grain we finally obtain 



L — ' m 



(16) 



The optical depth T\ of the dust is a function of k\, and 
k\ in turn is a function of T\ via the mass abundance 
Xi = Xi C^a)- In our models, the mass abundances Xi may 
change for two reasons. Firstly we are dealing with O- 
and/or C-stars. Secondly for the same type of star (either 
O- or C-stars) the opacities and the relative abundances of 
the grains can change with the optical depth (see Sect.QJ. 
Since k\ and T\ are each other interwoven, an iterative 
procedure is required . 

We also need the cross sections of the radiation-dust 
interactions. The cross section for a single dimension of 
the grains a is a 9i (a) = Tra 2 Q ex t (i), where Q ex t (?) are 
the extinction coefficients. The total cross section a gi can 
be obtained integrating the cross section for a single di- 
mension over a given distribution of the dimensions of the 
grains (see Sect. |S] for a discussion about the choice we 
have made for the distribution) . The absorption and scat- 
tering coefficients, Q a bs (i) and Q SC a («) respectively, can 
be calculated at any given wavelength with the aid of the 
Mie theory 2 . The parameters to be specified are: the op- 
tical constants of the material composing the dust grains, 
the dimension and shape of the grains that are usually 
approximated to a sphere. Basing on this, the absorption 
and scattering coefficients, Q a bs and Q sca , can be derived 
from the complex dielectric function, expressed by the real 
n (i) and t he imaginary part k (i) of the complex index o f 
refraction (|Bohren fc Huffmanlfl983l iDraine fe Leelll984j) . 

1 The radiative transfer code "DUSTY" refers the optical 
depth 7a at a wavelength of reference. Similarly to lSuhl |l999, 
2000j, I2002I) . we assume as reference wavelength A = 10/um. 
This implies that the term k\ in the optical depth equation has 
to be evaluated at this wavelength, once the dust composition 
is fixed. 

2 The Mie theory, developed in 1908 by Gustav Mie to un- 
derstand the colors generated by gold particles suspended in 
water, provides a formal solution for the interaction of spheri- 
cal and homoge neous small particle s with the electromagnetic 
radiation dBohren fc Huffmadll983h . 



The expansion velocity o f the m atter can be derived 
from the IVassiliadis fc Woodl l)l993|) relationship correlat- 
ing the expansion velocity to the pulsational period 



v exp [km s" 1 ] = -13.5 + 0.056P 



(17) 



with the condition that the velocity is higher than 3 km s 1 
llVassiliadis fc Woodl [l993) . The pulsation period P de- 
pends on the pulsation mode of the star. 
A simple relation between pe riod, mass and radius (eq. (4) 
in IVassiliadis fc Woodl Il993[) is obtained supposing that 
AGB variable stars pulsate into the fundamental mode: 



logP = -2.07- 



TT~ (18) 



where the stellar radius R and mass M are expressed in 
solar units. 

Another variable that appears in eq. I|12fl is the dust- 
to-gas ratio S. We need to r e late t his ratio to the stel- 
lar parameters. lHabing et al.1 <|l994^ studying the depen- 
dence of the gas outflow velocity v exp at large distances 
from the star, found that v exp depends on three pa- 
rameters, that is the stellar luminosity L, the mass-loss 
rate M, and the dust-to-gas ratio 8. The observational 
data indicate that v exp increases at increasing any one 
of these parameters as also confirmed by the models of 
lHabing et al.1 l|l994h who also suggested a plausible cor- 
relation between v exp and the dimensions a of the grains. 
Therefore, v exp = v exp (L,S, M,a) is a function of four 
variables. However, v exp is found to depend only weakly 
on the mass-loss rate M except when M is small and the 
same weak dependence hol ds good for the grain dimension 
a. See lHabing et al.l (^994) for all details. Since our aim is 
to include the effects of the dust shell around AGB stars 
on the spectra of SSPs, stars with a thick shell and high 
mass-loss rates give the dominant contribution to IR spec- 
trum of the SSP. Th erefore, we adopt here the relation of 
lHabing et al 1 lll994h for high mass-loss rates, and neglect 
the dependence of v exp on M and a. The expression for 
v exp in presence of high mass-loss rates is 



oc L u - 6 8> 



0.3 cO. 5 



(19) 



Followi ng iBressan et alJ lll998h . who made use of the re- 
sults bv lHabing et al.1 l)l994j) . we recast v exp as 



15 [k 



L 



m s 



("10" I. 



0.3 



( 6 



\8agb 



(20) 



and eventually invert it to estimate 6. The factor Sagb 
depends on the kind of star under consideration. For the 
O-s tars we use Sagb = 0.01 as in IBressan et al.1 lll998l) 
and lSuhl l|l999|) . whereas for the C- rich stars we adopt t he 
value Sagb = 0.0025 computed bv lBlanco et all l)l998|) . 
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3. Mass- loss rates from AGB stars 

The rate of mass- loss along the AGB is a key parameter for 
the evolution of the stars in this phase because it affects 
their lifetime and luminosity. 

Despite the great effort devoted to clarify the role of 
mass-loss on the evolution of AGB stars, the mechanism 
of mass-loss itself is not completely clear. The determina- 
tions of the mass-loss rate from IR and radio data, are 
affected by significant uncertainties, due to the poorly 
known distances, dust-to-gas ratio, and expansion velocity 
of the sources. 

There are two nearly direct observational evidence for 
mass-loss: the continuum IR emission in excess to what we 
expect from a star with the typical effective temperature 
of an AGB, and the molecular rotational lines detected 
in emission. The IR emission observed is characteristic of 
a dusty shell that absorbs the light of the star and then 
emits it in the IR and radio range. The molecular lines 
come from the gas surrounding the star: their widths and 
time variations prove that the gas is flowing away from 
the star. 

It is now widely accepted and supported by hydro- 
dynamical models that large amplitude pulsations are es- 
sential for accelerating the mass outflow from the stellar 
surface of AGB stars until the gas becomes cool enough 
that heavy elements can condense into dust grains. The 
dust grains in turn absorb and scatter the radiation trans- 
ferring by collisions energy and momentum from the stel- 
lar radiation field to the g as so that the flow velocity may 
exceed the escape velocity ( Gilman 1972). Mass-loss grows 
with time until the so-called super-wind regime sets in, 
which quickly turns the star into a planetary nebula by 
stripping away all the envelope and leaving a bare core 
th at evolves to high temperat ures. 

IVassiliadis fc Woodl l| 19931) represented the above sit- 
uation with suitable analytical fits in which M expo- 
nentially grows with the luminosity up to the forma- 
tion of the planetary nebu la. Following the formalism of 
IVassiliadis fc Woodl l|l99.'fl we have adopted the relations 
below to express the rate of mass- loss prior and during the 
super-wind regime. The transition occurs at a period of 
about 500 days. Furthermore, once the super-wind regime 
has started the stellar envelope gets so rich in dust that 
the star is no longer observable in the visible range of the 
spectrum. 

Prior to super-wind, the mass-loss rate is expressed by 



logM = -11.4 + 0.0123P 



(21) 



where M is in Mq yr 1 and the period is in days. During 
the super-wind regime M is assumed to be given by 



M = 



1 



L 



(22) 



which describes a wi nd driven by radiation pressure 
l|lvezic fe Elitzurlll994|) . 

As already pointed out bv IVassiliadis fe Woodl l)l993(l 
long period variables detected in the optical are found up 



to periods of 750 days (top upper part of the AGB where 
the stell ar mass is ~ 5M Pl ) . In or der to take this into 
account, IVassiliadis fc Woodl l|l993h slightly changed the 
expression for the mass-loss rate of eq. <|21[) according to 



logM = -11.4 + 0.0125 x 



100 



M 



2.5 



(23) 



for stars with masses M > 2.5A/©. Once again M is in 
M©yr _1 , and the period is in days. The equations (|2*l"1) . 
(|22|l and (|23|l completely describe our mass-loss prescrip- 
tions for AGB stars. 



4. Formation and properties of the dusty 
circumstellar envelopes 

O-rich stars. O-rich AGB stars of M spectral type show 
two typical features at 10/zm and 18/xto either in absorp- 
tion or in emission depending on the optical depth of 
the surrounding envelope. These features are usually at- 
tributed to stretching and bending modes of Si — O bonds 
and O — Si — O groups and clearly probe the existence of 
silicate grains in the shell of matter around the star. 

Among AGB stars, the OH/IR stars are generally 
thought to represent the final evolutionary stage of an O- 
rich object, just before it evolves quickly into a planetary 
nebula. Conversely, Mira variables of M spectral type are 
thought to correspond to early or advanced AGB evolu- 
tionary stages of oxygen-rich stars. An oxygen-rich star is 
characterized by the ratio C/O < 1. Because of the strong 
triple bond between O and C in the carbon monoxide, it 
is believed that all C will be blocked into CO molecules 
and no C is available to the formation of dust grain with 
other elemental species of low abundance. In contrast, the 
fraction of O not engaged in CO reacts with other el- 
ements such as Mg and Si and forms various types of 
compounds. Therefore, an O-rich star is characterized by 
the presence of O and compounds like MgO, silicates and 
H2O. Eventually, these molecules can bind together to 
produce grains of silicates. 

The newly formed dust grains leave the envelope of 
the AGB stars and disperse into the interstellar medium, 
where they can be strongly modified by chemical and 
physical processes, such as collisions, interaction with en- 
ergetic photons, a ccretion and destruction of t he mantle, 
and so forth (see iLi fc Mavo Greenberdl2002l for a re- 
cent review and referencing of the subject). Therefore the 
opacities of dust silicates in AGB stars are probably differ- 
ent from those derived from laboratory measurements of 
terrestrial or meteoritic material. Furthermore, the classi- 
cal opacities of silicates of the diffuse interstellar medium 
cannot be used to mod el the dusty envelopes of AGB stars 
jOssenkopf et al 

There are two final points to consider: i.e. the consis- 
tency of the opac ities in use with the Kram ers-Kronig dis- 
persion relations I Landau fc Lifshit3ll96ft see e.g. ) which 
imply that the real and imaginary parts of the complex 
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dielectr ic function e (li) = e i + it i are not independent, 
see also lBohren fe Huffman! l|l983|) . and the consistency of 
the theoretical results with the observational ones. In our 
case, the spectra obtained from radiative transfer models 
of AGB stars surrounded by dusty shells must agree with 
the IR obse rvations of the s e star s. For more details on this 
subject see lDraine fc 

As already pointed out by ISuhl (1999), the dust opac- 
ities are often adjusted in such a way that observations 
(spectra) are reproduced, however, without checking for 
the Kramers-Kronig dispersion re l ation at the same t ime 
{Volk fe Kwoklll988i lGriffinlll993t ISuh fe Jonedll997t) . In 
some cases ijOssenkopf et al.l Il992j) the physical consis- 
tency of the input opacity is secured, but no comparison 
between observational da ta and the oretical predictions is 
made. A rare exception is lSuhl l|l999|) who presents optical 
constants that satisfy both the Kramers-Kronig relation 
and properties of O-rich AGB s tars better than in previ- 
ous studies. In addition to this, ISuhl ( 1999) gives two es- 
timates for the silicate opacities, one for the warm grains 
(obtained reproducing the emission feature at 10/im of 
OH/IR stars) and one for cold grains (obtained repro- 
ducing the absorption feature at 10/j-to always of OH/IR 
stars). The variation of the silicate opacity with the ob- 
servational constraint to be reproduced can be understood 
by considering that AGB stars can be surrounded by thin 
or thick shells and that the grain temperature is likely to 
be cooler in thick shells than in thin shells. A simple ar- 
gument supporting this possibility is that in AGB stars 
surrounded by a thin shell of matter a large part of the 
dust will be at high temperatures, whereas in AGB stars 
with a thick circumstella r shell, dust will be at lower tem- 
peratures. Basing on this,|Suh ( 1999), argues that different 
kinds of grain (opacities) are required at varying optical 
depth of the shell: cold silicates are suited to thick shells 
with the lOfim feature in absorption, whereas warm sili- 
cates are more appropriate to thin shells with the lOfim 
feature in emission. The largest difference between the 
opacities of the two species occurs a t A > I3um (see the 
top panel of Fig.G) and according to lSuhl l(l999h the tran- 
sition between the two physical situations occurs at the 
optical depth 7io = 3 and A = lO^tm. The same values 
are adopted here. 

The mid-infrared spectra of oxygen-rich circumstellar 
envelopes are dominated by the silicate bands at about 10 
and 19/im and so much less attention has been paid to the 
other IR spectral regions until the advent of IRAS and ISO 
in particul ar. However, in ma ny IRAS LRS spectra of Mira 
variables ijSloan et al.lll996j) an emission band at about 
13fim has been detected tha t could be due to aluminum 
oxide l|Begemann et al.lll997h. but its o rigin is still a mat- 
ter of vivid debate (Posc h et alJ 119991) . Other weak sili- 
cate features around 10 and 18 \im have been detected by 
the ISO-SWS observations which can be attributed to the 
presence of oxide particles. In many AGB stars with high 
mass-loss rates, ISO high resolution observations revealed 
the presence of prominent bands of crystalline silicates, 



for instance enstatite (MqSiOz) and forsterit e {Mg^SiOij 
l|Waters et al.lll996t IWaters fe Molsterlll999j) . 

All these features spurred the interest toward the com- 
position of circumstellar dusty sh e lls and grain formation 
theories lQa.il fc Sedlm avrl ll 999h . Em3 <l2002l) presented 
a dust model for the envelopes of O-rich stars in which 
not only amorphous silicates but also crystalline silicates 
are included . The a dopted opacity functions are those of 

I. Taeger et all l|l998h (the bottom panel of Fig. QJ. The 
key parameter controlling an envelope composed by amor- 
phous silicates and crystalline silicates is the so-called 
crystallinity parameter a, that is the ratio between total 
mas s of c rystalline silicates to the total mass of silicates. 
iH3(2002), examining various levels of crystallinity, found 
that for the same value of a, crystalline silicate features 
are stronger for optical depths smaller than 5 and weaker 
for optical depths larger than 10 (always referred to the 
optical depth at 10/im). A high crystallinity is required to 
produce strong features in models with high optical depth. 

The intensity of the crystalline silicate features is likely 
to correlate with the mass-loss rate from AGB stars. As 
a matter of fact these features a re stronger in OH/IR 
stars with high mass-loss rates ijSvlvester et all Il999|) . 
Because the observations show that the same features 
are not present in stars with low mass-loss rates, one 
may perhaps infer that these stars are also defici ent of 
crystalline silicates in their envelopes |Suhll2002h . The 
proble m, however, needs deepe r analysis because ISuhl 
(200 21) andlKemoer et al.l l|200lj) find contrasting results. 
iKemner et al .1 (|200lh even including the presence of crys- 
talline silicates in the dusty envelopes, fail indeed to re- 
produce the f eatures in q uestion. 

Following ISuhl l)2002l) . we adopt here the crystallinity 
parameter a — 0.1 for stars with low mass- loss rates and 
moderately optically thick shells of matter (Tio < 15), 
whereas for OH/IR stars with high mass-loss rates and op- 
tically thick shells (Tio > 15) we prefer the value a = 0.2. 
It is worth recalling that a = 0.2 is fully sufficient to repro- 
duce the prominent crystalline features shown by OH/IR 
spectra. Finally, in all the models the relative contents 
of enstatite (MqSiOz) and forsterite {Mg2SiOi) are the 
same as m ISuhl l)2002|) . i.e. we adopt the same relative 
contribution for the two components. 

C-rich stars. C-stars are the evolutionary descen- 
dants of Mira variables. Through continuous dredge-up 
of carbon into the envelope during the thermally puls- 
ing AGB phase, these stars eventually reach a carbon 
abundance in the outer layers larger than that of oxygen 
(C/O > 1). When this occurs the formation of O-rich dust 
ceases and it is replaced by that with C-rich compounds: 
the C-star phase begins. Therein after, the continuous for- 
mation of C-rich dust makes the envelopes of these stars 
more and more optically thick. By loosing mass at very 
high rates, they get enshrouded by thick envelopes that 
absorb and scatter the UV-optical radiation into the IR 
and radio range. 

Almost all these stars show an emission feature at 

II. 3/im, due to silicon carbide (SiC), whose presence 
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Fig. 1. Upper panel: extinction coefficients for two crys- 
talline silicates i.e. forsterite from Jaeger (2002, privat e 
communication) and enstatite from iJaeeer et alJ (1998). 
Bottom panel: extinctio n coeffi cients for cool and warm 
silicates taken from Suh (2002). 



was predic ted by Gilmanl (Il969l) and observationally con- 
firmed by lHackwelll l|l972h . Many infrared observations 
have brought into evidence and confirmed the presence of 
several types of dust grains in C-rich AGB stars: amor- 
phous carbon (AMC), silicon carbide, and magnesium 
sulphide (MgS) to mention the three dominant types. 
ISiihl ( 1999, 2000) has derived new opacities for the AMC, 
that are consistent with Kramers-Kronig dispersion rela- 
tions, and has also reproduced the observational data by 
means of suitable models of the dusty shells around AGB 
stars (see Fig. |5J| . The models are based on the new opac- 
ities and a complete treatment of the radi ative transfer, 
thus much improving upon prev ious st udies IIBlanco et alJ 
Il998l IGroenewegen et al.l ll99Sft . The i|l 999L l2000h 

models are characterized by two components, SiC and 
AMC (for typical models of carbon stars, see for example 
lLorenz-Martins et al.l l)200l|) and references therein). The 
30^m MgS feature is observed in a wide range of sources, 
going fr om low mass- loss carbon stars to planetary nebu- 
lae l|Honv et alJl2002|) . Because it is observed in C stars of 
large luminosity it is expected to significantly contribute 
to the integrated spectrum in the 30/xm range. However, as 
the optical constants of MgS are n ot measured for a suffi- 
cient ly wide range of wavelengths llSuhll2000t iHonv et alJ 
2002) the presence of MgS is neglected here. 

The grains of silicon carbide can be separated in 
two groups, hexagonal o r rhombohedral aSiC and cu- 
bic 0SiC. Many authors l|GroenewegerJ ll995llSpeck et alJ 
Il997t iBlanco et alll998|l argue that in order to fit the IR 
spectra of C-stars, aSiC is required sometimes together 



with (3SiC. However, all studies of meteoritic SiC grains 
have revealed the presence of only 0SiC. Since it is ther- 
modyn amically un likely that aSiC transforms itself into 
(3 SiC ()Bernatowiczlll997| . there seems to be a clear dis- 
crepancy between the aSiC suggested by the AGB spectra 
and the 0SiC indicated by the meteoritic data. As pointed 
out by ISpeck et alJ (^)99) the discrepancy could be due 
to inadequate KBr corrections 3 to the laboratory spectra 
so that a plausible way out can be found. 

AMC and SiC affect the spectrum in a different fash- 
ion: while the effects of AMC propagate over the whole 
spectrum, those of SiC are li mited to the Hum feature 
as ob servat ionally indica t ed. lLorenz-Martins fc Lefevrd 
l|l994|) and IGroenewegen! l|l995|) suggest that the ratio 
SiC to AMC decreases at increasing optical depth of the 
dusty enve lope. This is also con firmed by the models of 
ISuhl l)l999j) . Therefore, following ISuhl l)l999l) prescription, 
we assume here that the chemical composition of the dust 
changes with the optical depth of circumstellar shell. For 
optically thin dust shells (Tig < 0.15), the strong 11/im 
feature requires about 20% of SiC dust grains to fit the 
observational data for AGB; for dust shells with interme- 
diate optical depth (0.15 < T w < 0.8) about 10% SiC 
dust grains are needed, whereas for shells with larger op- 
tical depths, in which the llfim feature is either much 
weaker or missing at all, no SiC is required. Following 
[ Suhl dl99£ ), wh o used the optical constants of aSiC by 
IPfeourieT i 1988|) to calculate the opacity for SiC and re- 
produce the spectra of C-stars similarly, we adopt here the 
same data (see Fig.|2J. 

5. Parameters of the shell model 

An important parameter of the dusty shell models is 
the distribution of the grain dimensions. Nowadays many 
recipes have been proposed to reproduce the properties 
of t he dust grains: the y go from the classical power- 
law (iMathis et alJ 1977 ) to a log-normal distribution of 
IWeingartner fc Draind pOOlfr . In this study we adopt the 
simple distribution expressed by the Dirac delta function 

n (a) — S (a — cto) 

with ao = 0.1/im, where a is the dimension of the grains. 
Despite its simplicity, this type of distribution is fully ad- 
equate to our purposes and it has already been widely 
used t o model the dusty shell s of AGB stars llHash motol 
19951 livezic fc Elitzurl Il995l Ivan der Veen et alJ Jl995[ 



Suh fc ,Toneslll997t ^119971119991 Ivan Loon et al.lll99 



3 iDorschner et all <ll97&t) have studied the effects of the a 
dispersive medium made of potassium bromide (KBr) on the 
spectra of small quantities of silicates with the grain size. They 
find a shift of the spectral features emitted by the grains using 
either a KBr-matrix or a KBr-pellet in which the small grains 
are embedded. They explain the shift as due to an effect caused 
by the KBr-matrix that needs to be correc ted before us i ng the 
spectra obtained with the KBr-pellet. See lSpeck et alJ (1 1999ft 
for more details about the problem of the KBr correction. 
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Fig. 2. Absorptio n coefficien ts for amorphous carbon 
(dashed li ne) fromlSuhl ll2000|) and silicon carbide (solid 
line) from lPegouria l|1988|) . 



Isnhll2nnd Ekv. 3). Support to this kind of approximation 
is bv lBagnulo et alJ 1)19951 ) who find that a delta- function 
fits the spectrum of IRC+10216 much better than a power 
law. 

Since different types of grains are present in the dust 
mixture of the circumstellar shell, this could imply dif- 
ferent values of the grain temperature at the same dis- 
tance from the central source. As a consequence of it, the 
condensation temperature at inner boundary of the shell 
could vary with the type of grain. In the case of shells 
of the C-rich stars, for which we use a mixture of AMC 
and SiC, Td is assumed equal to lOOOX for both types of 
grain: in reality AMC and SiC should have their own tem- 
perature and temperature profile across the shell. This of 
course would add additional parameters to the problem, 
i.e. two condensation temperatures and two dust-to-gas 
ratios, thus rendering the whole problem much more com- 
plicated and beyond t he technical capability of the cod e 
DUSTY we are using l(Nenkova. Ivezic. fc Elitzurlll999h . 
The code indeed treats a mixture of grains using a pseudo- 
grain whose optical properties are the average of those of 
individual species and of course deals with a single value 
of Td- Fortunately, in the case of C-rich stars, the abun- 
dance of SiC is very low so that the approximation to a 
single value of Td is physically acceptable. By the same 
token, we may also assume the same value of T d =100 K 
for this type of shell. As first noticed bv ISubl l)l999|) . if 
in an expanding stellar envelope dust grains condense in 
amorphous form and then get a crystalline structure on 
a short time scale, it is likely that two different types of 
grain have the same temperature distribution. 



Another important parameter is the density profile of 
the grains across the shell. Many different density laws 
have been suggested: the simple power-law p oc r~ x is a 
popular choice even if the exponent x varies from author 
to author. The law p oc r~ 2 is often adopted. Equally good 
alternatives are p oc r" 1 or p oc r~ 3 . The latter seems more 
suited to describe the case of O-stars with thick envelopes. 
Furthermore, models with radiation pressure or models 
with pulsations and shocks predict p oc r -2 . Finally, a 
region of finite thickness with enhanced density is added 
to the distribution p oc r~ 2 to somehow describe stars 
that are in the super- wind phase l|Suh k. Joneslll997t ISuhl 
Il997j) . In this study we simply a dopt the power-law p oc 
r" 2 as in the series of papers bvlSuhl (fl999l l200d l2002h 
from which we also take the opacities of the dust grains. 

However, using different power-laws for the distribu- 
tion of the matter in the shell implies that the same over- 
all optical depth is reached at different radii. In order to 
evaluate the uncertainty on the optical depth T\ caused by 
using different density profiles across the shell. We start 
from eq. (@J, insert the relation pd(r) = 5p(r) where 8 
is a constant, suppose for the sake of simplicity that the 
product k\(r)S does not depend on r, and finally insert 
the generic radial dependence for the density po ( r / r o) " 
(where po and ro are suitable scale factors). The optical 
depth down to a generic radius r is 



T x (r) = (k x (r)5) / p{r)dr 

= (k x (r)6)p r% [ r-dr 



(24) 



At any given radius r the variation caused by using dif- 
ferent power-laws is 



A 



- n dr 



(25) 



Let us indicate with T\ , I2 and I3 the three integrals for 
three typical power-law exponents n = 1, 2 and 3, respec- 
tively, pose r = ari n , and finally assume the scale length 
ro to be the inner radius of the shell r„. We immediately 
get 



li = r m lna 1 2 



1 



1 



J 3 ~ 0.5r r , 



1 



1 



The factors multiplying Ti n give an idea of the differences 
due to the power-law exponent n. The net consequence 
of different mass density profiles is that the more matter 
is stored at large radii and hence lower t emperatures the 
flatte r is the spectrum and vice- versa fsee llvezic fc Elitzurl 
119971 for more details). 

6. Uncertainties of the shell model 

In this section, we try to discuss the uncertainties arising 
from the two main components of the adopted model for 
the dust shell. 
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Optical depth. Differentiating cq. (|12fl we get 



AS AM 

S M 



Afc> 



An., 



(26) 



In the following we try to estimate the various terms 
composing this relation. If we look at eq. 1)20 [). this one 
well represents the link between luminosity of the star 
Lj., outflow velocity v exp and dust-to-gas ratio 5, because 
ijHabing et alJl9 94) if L« or S increase, also v exp increases, 
but the exponents that rule the relation of proportion- 
ality are uncertain. W e adopt the depe ndence proposed 
bv lHabing et al 1 (ll994h . but for example Ivan Loonl <|2000|) 
proposed a slightly different choice for the L dependence. 
If we write eq. ()2U|I as v exp oc L"5@ we get 



AS_ 
T 



1 Aw. 



13 v, 



cxp 



cxp 



a AL 
+ J3~L 



(27) 



where a and (3 are 2 and 0.6 in our model. The errors 
Av exp and AL are small, because equations (|T7|) and (jT%|l 
are sufficiently good and so, if we consider a large spread 
of a and (3 and typical values for the expansion velocities 
and luminosities of AGB stars, we get AS/ 8 ~ 1 in the 
worse cases. 

The uncertainty on the rate of mass-loss can be es- 
timated from eq. (|2T)l prior to super wind and eq. 
during the super-wind. Differentiating eq. I|21|l we obtain 



AM 
M 



1 



logioe 



0.0123AP ~ 0.034AP 



(28) 



which for the error on the period coming from (|18fl yields 
that AM / M is negligible. The uncertainty on eq. lll'L'l) is 
simply 



AM 
M 



AL 
~L 



Av. 



cxp 



u cxp 



(29) 



which is also negligible. 

Estimating the term Ak\/k\ is a cumbersome affair be- 
cause k\ depends on the mass abundances \i of the grains, 
which in turn are functions of the optical depth 7~\ and of 
the cross sections a 9i s of the dust-radiation interactions, 
which depend in its turn on the grain dimensions a and 
on the extinction coefficients Q e xt(i)- No simple way of 
checking the uncertainty can be found, but referring to 
the original sources of data for extinction properties. 

As already mentioned, the uncertainty on the expan- 
sion velocity v exp is simply estimated from eq. (|17H , which 
yield a negligible Av exp /v exp . 

Finally, the uncertainty on the inner radius of the shell 
can be derived from differentiating eq. (|14f> 



= 2- 



AT d 



AL 
> 

' L 



(30) 



For Td we have many different values in literature: vary- 
ing it in the realistic range from 1000 to 1500if and 
taking the typical luminosity for AGB stars we estimate 



At i n I ri n ~ 1. We can finally conclude that the terms 
Ar i n I Ti n and AS/S are the main source of error in the 
calculation of the optical depth T\. For individual stars 
the error may be as large as AT\/T\ ~ 2 which is accept- 
able for the purposes of our study. It will be smaller as 
soon as better determinations of the relationship v exp (L*, 
S) will be available. 

Chemistry. Although the list compounds we have con- 
sidered in modelling the chemical composition of the dust 
shells is certainly incomplete, an effort has been made to 
include as many molecules as possible. Many other com- 
pounds are of course possible and should be included in 
the list, but we need deeper studies able to discriminate 
for which AGB stars these compounds have to be in- 
cluded. Ho wever, compar e d to s imilar studies on the same 
topic, e.g. iBressan et alJ l)l998l) . our treatment is a sig- 
nificant step forward. To conclude we are confident that 
our description of the shell composition is adequate to the 
present aims. No doubt that this aspect can be improved. 

Final remarks. Even if the uncertainty is as large as 
AT\/T\ ~ 2, the situation is, however, not as bad as it 
may seem, because the largest effect of the dust shell on 
re-processing the radiation coming from the central star 
underneath is for the optically thick case, in which the 
exact value of the optical depth is less of a problem. 

7. Theoretical SSPs: basic ingredients 

The monochromatic flux of an SSP of age t and metallicity 
Z at the wavelength A is defined as 



rM v {t) 

ssp x (t,Z)= / f x (M,t,Z)$(M)dM 

J Ml 



(31) 



where f\ (M, t, Z) is the monochromatic flux emitted by 
a star of mass M, age t, and metallicity Z; <f> (M) is the 
initial mass function (IMF); Ml is the mass of the lowest 
mass star in the SSP whereas Mu(t) is mass of the highest 
mass star still alive in th e SSP of age t. For the IMF we 
adopt the ISalpeterl l|l955|) law expressed as 



dN 
dM 



$(M) = AM~ 



where x=2.35 and A is a normalization constant to be 
fixed by a suitable condition (SSPs for other choices of 
the IMF can be easily calculated). 

Because the flux ssp\ (t, Z) is calculated by integrat- 
ing equation l|31|) along an isochrone, one has to know the 
luminosity, T e ff and gravity of the stars of mass M, age 
t and metallicity Z lying on the isochrone. For any star 
of mass M, age t and metallicity Z, the relationship be- 
tween luminosity, T e ff, gravity with the age t is derived 
from a library of stellar models, the flux f\ (M, t, Z) from 
a library of stellar spectra. Finally one has to adopt a pre- 
scription for the amount of mass lost by a star of mass M 
in the course of its evolution: this could occur all over the 
evolutionary history such as in case of massive stars, or 
at the end of the RGB passing from the tip of the RGB 
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to the HB or clump as appropriate, and during the AGB 
phase as in the case of low and intermediate mass stars. 
All detai ls of the isoch r one co nstruction technique can be 
found in lBertelli et alJ l|l994h . 



7.1. Libraries of stellar models and stellar spectra 



In our st udy we adopt the isochroncs 
by iTantalol l)l998j) (anticipated in the 
data base for galaxy evolution models by 
iLeitherer. Alloin. von Alvensleben. Gallager. Huchra. fc et 
(1996)). The underlying stellar models are those of the 
Padova Library. T hey are shortly referred to as the 
iBertelli et alJ l)l994() stellar models (see references therein 
for more details). These stellar models are calculated 
wit h convective oversh ooting and are amply described 
by IBertelli et alJ l(l994h so that no detail is given here. 
In the present study only three chemical compositions 
are considered, namely [Y=0.240, Z=0.004], [Y=0.250, 
Z=0.008], [Y=0.280, Z=0.02], for which the transition 
luminosities for AGB stars pas sing from the O-ric h to 
the C-rich regime calculated by iMarigo et al.1 l|l999h are 
available (see below). 

In the isochrones, all evolutionary phases from ZAMS 
to the end of the TP-AGB or C-ignition are included, 
as appropriate to the mass of the last living star of age 
t. The rate of ma ss-loss during t he RGB stages of low 
mass stars is from iReimerd l|1975|) with r\ = 0.45. During 
the AGB phase the rate of mass-loss is according to 
IVassiliadis fc Woodl I 19931) . as already described in Sect. 01 
The initial masses of the stellar models go from 0.15 to 100 
Mq and the ages of the isochrones go from 3 Myr to 20 
Gyr. 

The library of stellar spectra is from iLeieune et alJ 
which stands on the Kurucz (1995) release of 
theoretical spectra, however with several important im- 
plementations. For T P ff < 3500 °K th e spectra of 
dwarf stars by lAllard k. Hauschildtl ( 19951) are included 
and for giant stars the spec tra by Fluks et al.l l)l994h 



-6 



-5 



Z = 0.004 



al. 



and iBessell et al.1 Jl989L Il99lh are considered. Following 
iBressan et alJ l|l994l) . for T eff > 50000 °K, the library 
has been extended using black body spectra. 

For ages younger than about 100 Myr the SSP s we are 
going to present will be equal to the old ones of ITantalol 
l)l998|) . because no AGB stars can be formed, whereas 
for older ages the effect of the circumstellar envelopes in 
processing the radiation coming from the central star is 
taken into account. 



7.2. Transition luminosities of the AGB phases 

Owing to the different composition of the ejecta much 
affecting the chemistry and physics of the circumstellar 
shell, the transition from the oxygen-rich M-type stars 
(C/O < 1) to the carbon-rich C-type objects (C/O > 1) 
deserves particular care. 




M=l.lM r 



Fig. 3. The range of bolometric magnitudes spanned by 
AGB stars with Z=0.004 and the transition bolometric 
magnitudes from M- to C-star as function of the initial 
mass. The dashed lines show the start of the AGB phase. 
The solid lines show the end of it. The shaded area is the 
luminosity interval in which C-stars are fo rmed. The thick 
lines (dashed and so lid) are the models bvlTantalol (^998) 
calculated with the IVassiliadis fc Woodl l)l993h prescrin- 
tion for mass-loss during the AGB phase. The thin lines 
(dashed and s olid) s how the same but for the models by 
Marig o et al. (1999). The region of C-stars is according 
to the IMarigo et al (^999) models. The two vertical lines 
show the minimum mass for the formation of C-stars and 
the maximum mass M up iov the occurrence of the AGB 
phase. 



In view of the discussion below, it is useful to sum- 
marize the evolution of AGB stars eventually becoming 
C-stars. The situation is illustrated in Figs. E3 El and 
for three different values of the metallicity, namely 
Z=0.004, Z=0.008 and Z=0.02 respectively. Each plot 
shows the bolometric magnitude of TP-AGB stars as a 
function of their initial mass Mi from the beginning to 
the end of the AGB phase when the outer envelope is 
completely removed by mass-loss. The thick dashed and 
solid line bind the luminosit y attained by the ste llar mod- 
els of th e Pad ova Library ijBertelli et al.l fl994) used by 
ITantalol J 1998) assum i ng th e prescription for mass-loss 
of lVassiliadis fc Woodl dl993h . The thin dashed and solid 
lines show the sam e according to the m ore recent models of 
TP-AGB stars bv lMarigo etaD (|l999|) . The shaded areas 
co rrespond to th e lumi nosity interval in which, according 
to IMarigo et alJ l|l999ft . C-stars develop. Exam i ning the 
three diagrams and the models bv lMarigo et al 1 lll999h m 
detail we note that: 
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Fig. 4. The same as in Fig. 01 but for the metallicity Z = 
0.008. 
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M, 

Fig. 5. The same as in Fig. but for the metallicity 
Z=0.02. 



(i) Over the mass interval in which the AGB phase de- 
velops, the formation of C-stars does not occur for initial 
masses lower than a certain limit indicated by the verti- 
cal line (the limit mass increases at increasing metallicity). 
This is simply caused by the very early loss of the envelope 
before any change in the chemical composition of the outer 
layers may take place. Therefore, in old stars only O-rich 
shells of dust with C/O < 1 are possible. At larger initial 



masses, C-stars and C-rich dusty shells do occur over a 
wide range, which may extend up the so-called M up , i.e 
the maximum value for the occurrence of the AGB phase, 
or slightly below this. In the range of the massive AGB, 
the occurrence of envelope burning may indeed delay the 
formation of a C-star, and in some cases may turn a Ci- 
st ar into a late O-rich star. More precisely, in this mass 
range a C-star is comprised between two stages, in which 
an AGB star appears as O-rich. The late O-rich phase is 
favored at low metallicities and disappears for Z > Zq. 
(ii) As a consequence of the above trend, for the majority 
of C-stars (say up to Mj < 3.5M Q ) the maximum luminos- 
ity just coincides with the end of the AGB phase, whereas 
for initial masses larger than Mj > 3.5Mq the maximum 
luminosity of C-stars may be fainter than the maximum 
luminosity of the AGB phase. 

Giv en these premises, w e need to incorporate the re- 
sults bv lMarieo et al.l l)l999(l . which provi de the transitio n 
luminosity into the isochrones (SSPs) bv iTantalol (^98). 
The comparison of the initial, transition, and maximum 
luminosities in the two sets of models reveals that: 

(a) For Z = 0.004 the luminosity at which the TP- 
AG B phas e begins is nearly th e same for both ITantalol 
(199$ and lMarieo et al.l (fl9 99). whereas the termination 
luminosity is lower in ITantalol l|l998| ) . Finally C-stars are 
formed for initial masses larger than about 1.1M Q . Owing 
to the lower termination luminosity, the proportion of C- 
sta rs with respect to that of the O-rich ones predicted by 
the ITantalol |l998) will be somewhat smaller compared to 
the value expected from lMarigo et all l|l999l) . 

(b) Similar considerations apply to case with Z = 0.008, 
the only difference being that no C-stars can be formed for 
initial masses lower than about 1.2M . We note that the 
ra nge of b o lomet ric magnitude predicted b y the models 
of iTantal o (1998) is narrower than that bv iMarigo et alJ 
( 199^), with lower maximum and higher initial luminosi- 
ties for the AGB phase. This reduces the allowed range 
for both M- and C-stars. Nevertheless, to a first ap- 
proximation the percentage of C-stars with respect to 
that of O-rich stars remains similar to that predicted by 
IMarigo et ail ljl99flh . 

(c) For Z=Zq, th e initial and ma xim um luminosity of the 
TP- AGB phase in lTantalol ljl99Sh and lMarigo et alJ l|l999h 
are nearly t he same s o tha t th ere would be no diff erence 
using either ITantalol l)l99Sj) or IMarigo etaH l|l999|) mod- 
els. Now the formation of C-stars occurs for initial masses 
larger than 1.5M Q . 

(d) Finally, for the sake of simplicity the possibility of a 
late O-rich phase is taken into account only for Z=0.004 
and Z=0.008 and is ignored at higher metallicities owing 
to the very small difference between the maximum lumi- 
nosity of the AGB phase and of the C-star regime. See for 
instance the case of Z=0.02 in Fig. 

Basi n g on the above considerations the use of the 
Tantalol l|l998f) isochrones is fully adequate to our aim 
of studying the effect of the dusty envelope on processing 
the radiation co ming from the c entral object. Therefore 
we will adopt the ITantalol l|l998|) models as source of our 
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SSPs and use the models bv lMarigo et al, | lll999h only to 
fix the luminosity at which the transition from O-rich to 
C-stars o ccurs . Work is in progress to replace the SSPs by 
iTantalol {l998) models with new ones fully incorporat ing 
the TP- AGB stars as modeled bv lMarigo eTaD l|l999h . 

7.3. A remark on the libraries of stellar spectra and the 
evolutionary models of AGB stars 

The general accuracy of the final result depends on the 
adopted libraries of stellar spectra, in particular the cor- 
relation between the theoretical parameters T" e //, gravity 
and chemical abundance and the associated SED, and im- 
portant details of the evolutionary models for AGB stars, 
in particular the dredge-up episodes altering the surface 
chemical abundances, the efficiency of mass- loss terminat- 
ing the AGB evolution, the path of AGB stars in the HR- 
Diagram eventually determining the T e ff and luminosity 
of the AGB stars and the associated SEDs in turn, and fi- 
nally the accuracy of the isochrone construction and of the 
stellar models underneath. To summarize and review in 
detail all the sources of the various items above would lead 
to too a lengthy discussion which goes perhaps beyond the 
scope of this study. Suffice it to mention, that the libraries 
of stellar spectra, isochrones, SSPs in usage here have been 
successfully tested against, somehow tailored to match, 
and applied to study the broad band UBVRI photometric 
data of star cluste r s (ste ll ar popul a tions i n general), e.g. 
lGha,rdi Bertelil JToijft. lOira.rdil Jl99l. ICa.rra.ro et al] 
l)l999(l . iTantalo et al.l l|1998() . As far as the evolutionary 
models of AGB stars are concerned, the same considera- 
tions hold good, because they have been calibrated on the 
observational data, e.g. t he luminosity funct ion of C-stars 
in the Magellanic Cloud l|Marigo et al.ll 999). To conclude 
the libraries of stellar spectra and evolutionary models of 
AGB stars wc have adopted arc fully adequate to our pur- 
poses. It is, however, long known that the IR colors of the 
classical models of AGB stars disagree with the observa- 
tional data. In addition to improving the physics of the 
models, for instance by including the effect of C-ri ch com- 
positio n on the opacity of the outermost layers, see lMarigol 
( 2002) which would yield AGB models of much cooler T e ff 
during the C-rich phase (these models are not yet incor- 
porated in this study), the dust shells are expected to 
strongly affect the IR colors and to bring them to better 
agree with the observations. This indeed is the aim of this 
study. 

8. Infrared spectra of SSPs 

In Fig.0we plot the SED of dusty SSPs for different values 
of the age limited to the case with metallicity Z=0.008 and 
in Fig. we show the same but for the classical SSPs in 
which the effect of dust is neglected. The displayed SEDs 
are for ages in the range 0.25 to 10 Gyr. A better view of 
the difference brought by dust is shown in Fig. [SJ in which 
the new and old SEDs are compared for a few selected 
ages. The differences are remarkable. First of all, in the 




1 10 100 

A (/xm) 



Fig. 6. Integrated SEDs F v vs. A for the SSPs with Z = 
0.008, ages from 0.25 to 10 Gyr and the inclusion of dusty 
circumstellar envelopes in AGB stars. From the bottom to 
the top the displayed ages are : 10, 7.5, 5, 4, 3, 2, 1.5, 1, 
0.8, 0.6, 0.5, 0.4, 0.35, 0.3, and 0.25 Gyr. 

old SSPs the spectra do not extend into the medium and 
far IR (MIR and FIR, respectively), but sharply decline for 
wavelengths longer than about 3 — 4 pm. In contrast, the 
spectra of the new SSPs extend toward long wavelengths 
and the flux is considerable also in the MIR and FIR. 

The differences start at about lpm and in the IR range 
up to 3 — 4 pm the flux of dusty SSPs is lower than the 
old one: this is a consequence of the fact that dusty en- 
velopes shift the emission of M and C stars toward longer 
wavelengths. The amount of energy shifted to longer wave- 
lengths is larger for the young ages, because more massive 
and luminous AGB stars are present. 
It is worth noticing the different IR spectrum of the new 
SSPs, the evolution of the features at 11.3 pm and 9.7 pm 
in particular. For young ages, e.g. the SSPs of 0.3 Gyr 
and 1 Gyr, the spectrum does not exhibit features due 
to crystalline silicates, because the C-stars dominate (the 
11.3/£m feature of SiC is indeed prominent); for interme- 
diate ages, such as 3 Gyr, the O-stars of low optical depth 
influence the spectrum and the 9.7 pm feature can be seen 
in emission. For older ages, from 5 Gyr onward, the O-stars 
dominate, so the spectrum becomes more articulated and 
the features due to crystalline silicates start to appear at 
long wavelengths in the IR. 

In Fig. we plot the detailed evolution of the SiC 
and Si — O stretching mode features at increasing age. 
At young ages (0.3, 0.5 Gyr) there is only the feature at 
11.3pm of the SiC; at about at 1 Gyr, the 9.7pm of 
Si — O starts to appear; in the age range 1 to 2 Gyr the 
two features overlap; finally, at older ages the 11.3 pm of 
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0.3Gyr 




Fig. 8. Detailed comparison of SEDs F\ vs. A for the old 
(dashed lines) and new SSPs (continuous lines) with Z = 
0.008. Only four ages are considered as indicated. 



feature is slightly in absorption rather than in emission as 
at older ages. The explanation of it can be seen in Fig. [5J 
For the youngest ages, the luminosity interval in which 
C-stars appear becomes very thin, and the stars at the 
AGB tip, with the highest mass-loss rate and the highest 
optical depth, are O-stars. So the SSP spectrum becomes 
dominated by the 9.7/im feature in absorption and not in 
emission, as expected in envelopes with high optical depth. 
A similar situation would also occur at lower metallicities 
(e.g. Z = 0.008), but only at ages younger than 0.2 Gyr 
as shown by the luminosity intervals for C-stars displayed 
in Fig. |H 

ISilva et afl l|l998h calculated spectra of SSPs with 
AGB dusty shells for chemical mixtures containing either 
silicate or graphite grains. They do not consider the more 
realistic situation in which both types of grain are present 
and do not include the transition of O-stars to C-stars 
(followed in some circumstances by a late transition from 
C-stars back to O-stars for AGB stars of very high mass). 
Thanks to the inclu sion of the transition luminosities of 
iMarieo et alJ l)l999|) . it is now possible to follow the IR- 
SED at changing the proportions of C-stars with respect 
to O-stars, which is expected to vary with the evolution 
of the AGB star population. 



SiC disappears, and only the feature at 9.7 \im of Si — O 
occurs. 

In order to illustrate the effect of the metallicity, we 
plot in Fig. El the SEDs of the case with Z=0.02 for the 
same ages of the case with Z=0.008 shown in Fig.[f)] The 
evolution of the spectrum is similar but we can observe an 
interesting change for the youngest ages where the 9.7 \xm 



9. Broad-band colors in the infrared 

To test how the presence of dusty shells around AGB stars 
would affect the SEDs and colors of individual objects, we 
derive IRAS and Johnson broad-band colors and compare 
them with the observational data for a selected sample of 
AGB stars of different type. In addition to this we calcu- 
late the integrated colors of SSPs in the Johnson system 
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Fig. 10. Evolution of the 11.3 /Ltm feature of SiC (dashed 
vertical lines) and the 9.7 fmi feature of stretching vibra- 
tions of Si — O bonds (dotted vertical lines) caused by the 
stretching modes of this molecule in SSPs of metallicity 
Z = 0.008 and ages going from 0.3 Gyr to 5 Gyr. 

and compare them with a sample of star clusters of the 
Magellanic Clouds. 

IRAS Colors. To derive the theoretical IRAS 
monochromatic fluxes we have convolved the SEDs of 
our SSPs with IRAS transmission curves, derived from 
the IRAS Explanatory Manual Supplement (on line 
at http:/ /space, qsfc.nas a. gov/ astro/iras/ docs/ exp.sup/) . 
Following iBediinl 
fined as 



the monochromatic flux is de- 



F A = 



/ F A $ A dA 



/(f) 



(32) 



where $a is the instrumental profile, F\ is the theoretical 
flux distribution, F£ is the assumed flux distribution, that 
for IRAS bands is oc A -1 , and Fq A is the assumed flux 
distribution referred to the central wa velength of the band 
as defined in Neugebauer et all l|l984|) . If the source has a 
profile with a shape different from the dependence oc A -1 , 
one has to correct the flux using the so-called K factor 
given by 



K = 



j (^) $udv 



(33) 



and to divide the monochromatic flux resulting from eq. 
<|32l) by this factor. Finally, the IRAS colors are given by 



[12-25] = -2.5 log 



F(12) c 



2.5 log 



f (25) 
F(25) 



(34) 



[25 - 60] = -2.5 log 



F(25) 



2.5 log 



F(60) 
F~WT 



(35) 



where F(12), F (25), and F (60) are the IRAS theoret- 
ical fluxes in Jansky (Jy), and the constants of calibra- 
tion F (12) =28.3 Jy , F (25) Q =6.73 Jy, and F (60) =1.19 
Jy are taken from the IRAS PSC Explanatory Manual 
Supplement (1988) 4 . 

Our template sample of O-rich AGB stars (Miras, 
Semi- Regular, Long Period Variables, OH/IR etc, indi- 
cated by the filled circles) and C-stars (open circles) 
is shown in the IRAS two color plane [25-60] vs [12- 
25] of Fig. The location of the stars in this di- 
agram is widely used to establish t he nature of the 
sources. Ivan der Veen fc Habinel 1 19881) have defined a re- 
gion named the VH-window occupied by late-type stars. 
In addition to this, many other diagnostic plots are used 
to discriminate between O stars and C stars. The two- 
colors plots are constructed combining in various ways 
IRAS fluxes w i th oth er IR pass-bands, see for instance 
lEpchtein et alJ l)l987[) . In this diagram we also plot the 
colors of the AGB stars. In order to show the whole range 
of colors spanned by these stars we make use of the SSPs: 
more precisely the colors displayed by an SSP would cor- 
respond to AGB stars of the same age but different initial 
mass. By varying the age of the SSPs we may cover the 
whole range of colors of AGB stars of any initial mass and 
age. 

In Fig^] models of classical AGB stars, i.e. with no 
dusty shells around, would span a small range of colors ap- 
proximately centered at about around [12 — 25] ^ and 
[25 — 60] i2 (the encircled area visualizes the region cov- 
ered by two SS Ps with age o f 0.25 and 5 Gyr). As already 
pointed out bv ISilva et alJ l)l998j) . this happens because 
the SED of AGB stars, whose T e ff falls in the range 2500 
to 3500 K, is much similar to the Rayleigh- Jeans distribu- 
tion of black-body. The net result is that classical models 
of AGB stars fail to reproduce the observational distribu- 
tion of O and C-stars in the VH-window. 

Passing now to the new models of AGB stars with the 
dusty shell around, the situation is much improved. This 
is shown in Fig.^J where now the colors of the AGB stars 
(represented by two SSPs of given age) stretch across the 
whole VH-window. As we expect, the path in the two color 
plot of massive AGB stars (the young SSP of 0.25 Gyr 
shown by the dashed line) and low-mass AGB stars (the 
old SSP of 5 Gyr shown by the solid line) is different. The 
low mass AGB stars overlap only the region occupied by 
O-stars, whereas the more massive AGB stars jump into 
the region crowed by C-stars when the transition from 
O- to C-star occurs. This stage is marked by the dotted 
portion of the dashed line. 

However, it must be noticed that many stars in our 
sample have [25 — 60] colors that are significantly redder 
than predicted by the theory. 

4 We have also calculated the flux F (100) adopting the cal- 
ibration constant F (100) =0.43 Jy. No use of this quantity is 
made here. 



16 



L.Piovan, R. Tantalo & C. Chiosi: Shells of dust around AGB stars 




Fig. 11. The IRAS two color diagram [25 — 60] vs 
[15 — 25]. The O-stars have been sa mpled from dif- 
ferent sou rces: iLewis et al l l|l990ft , |Le_Scjuereii i ^^ll 
1 1992ft . iDavid et alJ il993ft . Blommaert et alJ 



lll993ft. IChengalur et all lll993ft . ILoup et all 1 1993ft . 
Xiong et alJ ill 994ft. I^nine et alJ ill 99,5ft . iLewisl ill 997ft . 

van Loon et all Jl998ft . Th e sample of C- st ars is 
from lEpchtein et all lll99f]ft. lEgan fc Leunel Jl991 
IChanl lll993ft. IGuglielmo et all lll993ft. IVolk et all 
lll992ft. iGroenewegenl (|l995ft . iGuglielmo et al 
IGuglielmo et all (|l997ft . The fluxes of the stars in the 
samples are obtained using the IRAS Point Source 
Catalogue. Two groups of coeval AGB stars of different 
mass represented by the SSPs with age of 0.25 Gyr 
(dashed line) and 5 Gyr (solid line) are compared to the 
data. The large circle marks the color range spanned by 
these AGB stars. Classical models of AGB stars are not 
able to fit the observational color distribution. 

This fact could b e easily explained by the so-called 
cirrus contamination l|lvezic &: Elitzurl ll995). In short the 
cirrus emission C\ may affect the IRAS fluxes because, ow- 
ing to the point-like nature of the sources, sky-subtraction 
may not be accurate enough and some contamination by 
the cirrus light can be present and significant, so that 
some additional correction of the data is required The 
cirrus emission Cioo can be estimated from IRAS quan- 
tity cirr3, i.e. the surface brightn ess at lOOum around th e 
point source: Cioo = 1.2-cirr3 Jw. llvezic fc Elitzurl l|l995ft . 
have shown that long- wavelength IRAS point source fluxes 
are unreliable when cirrS > (1 — 5) --F60- Basing on this, it 
is possible to clean the sample by removing all the sources 
with > where 1 < a < 5. In Fig.El are shown the 
same data of Fig. ll2l but in which all stars with c 'J r3 > 3 
have been removed. Now the theoretical colors much bet- 
ter reproduce the observed ones. In general at increasing 



Fig. 12. The same as in Fig. II II but in which two groups 
of new AGB stars of the same age and different initial 
mass are displayed. Like in Fig. ^2 these AGB stars are 
simulated by two SSPs with age of 0.25 Gyr (dashed line) 
and 5 Gyr (solid line). Along the line for the young age 
we mark with a dotted line the rapid transition from O - 
to C-stars. 



5Gyr 
0.25Gyr 




1 12-25] 



Fig. 13. The same as in Figs. Illlandll2lbut in which the 
data have been selected imposing the condition C p^ 3 < 3. 
All the stars with ratio £Hz2 greater than 3 have been 

^60 ° 

discarded because too much contaminated by the cirrus 
light. 



L.Piovan, R. Tantalo & C. Chiosi: Shells of dust around AGB stars 



17 



the cirrus correction, the agreement between theory and 
data improves, but too few stars are left over and the com- 
parison loses statistical significance. 

Another point to note is that our models of AGB stars 
of different mass do not actually cover the whole color 
ranges of the data even considering the cirrus correction 
above. It is likely that the parameters adopted in the ra- 
diative transfer problem are not fully adequate. There are 
several causes worth being examined in some detail: 

(i) Part of the disagreement can be due to the opti- 
cal depth of the envelope and to its dependence on the 
chemical composition. Recalling the definition of optical 
depth, this could vary from model to model because of dif- 
ferent compositions of the dusty envelopes. In this study 
we adopt different compositions of the dust for different 
intervals of optical depth, however we have not considered 
the possibility of different mixtures of silicate and carbon 
grains in envelopes with the same optical depth. T his is 
a point of weakness because llvezic fe Elitzurl l|l995l) have 
demonstrated that mixed compositions could explain the 
spread of the data (see for example their Fig. 6). 

(ii) Another possibility, is that the density profile of the 
matter is more complicated than a continuous power-law, 
as a consequence of discontinuous episodes of mass-loss of 
different intensity. Therefore the approximation p oc r~ 2 
could be too simpl e, actua lly only mi r roring average prop- 
erties. Em3 <|l997h and ISilva et all l|l998li have indeed 
shown that a narrow r egion of enhanced density, c aused 
by a super- wind phase IjSteffen fc Schonberne r 2000) trav- 
elling across the envelope, could be the cause of a large 
spread in the two color diagram. We plan to examine this 
problem in a forthcoming study. 

NIR Colors. Wc also compare our theoretical colors 
in near IR, such as [J — H], [H — K] and [K — L] with 
the data in the same pass-bands for a sample of O- and C- 
stars. The comparison is shown in Figs. 1141 and 1151 Three 
groups of AGB stars are shown: massive objects (the SSP 
of 0.2 Gyr, dotted dashed line), intermediate mass AGB 
stars (the SSP of 5 Gyr, solid line), and low-mass AGB 
stars (the 9.5 Gyr SSP, dashed line). While the classical 
models for AGB stars are confined in a small region and 
thus fail to match the data, the new models with the dusty 
shells stretch across the whole region crowded by the ob- 
servational data. Both in the [J — H],[H — K] and in the 
[J — K], [K — L] two color diagrams the fit is good, even 
if we still ha ve some problem t o fit the data of O stars. 
As noticed by Sil va et alJ l|1998h . problems in the fit could 
be due to the atmospheric models adopted for the M type 
spectra. 

Star Clusters. It might be worth of interest to 
compare the integrated broad-band colors of the SSPs 
whose AGB stars are enshrouded in the dust shell to 
those of a sample of star clusters. To this aim we have 
looked at the young globular clusters of the Large and 
Small Magellanic Clouds (LM C and SMC, respectively). 
IPietrzvnski fe Udalskil l|200Cf) have presented age determi- 
nations for about 600 star clusters belonging to the cen- 
tral part of the LMC. These clusters are younger than 1.2 
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Fig. 14. The two color diagram [H — K] vs [J — H] 
derived from the new models of AGB stars. The colors of 
the O H/ IR and Mira st ars ar e t aken from ILepine et alJ 
lll995l). IXionget al.1 jl994h . lOlivier et all ll200lh . 
IWhitelock fe Feastl <ll99HT whereas those of the C-stars 
are derived from lEpchtein et all lll990h. IGuglielmo e" 



lll99fllOuglielrno etal] ill 997ft. IGuglielmo et~aH <jic 
lOlivier et alJ l(200lh . All the data are properly cor- 
rected for extinctio n. In particular, for the database of 
ILepine et all l|l995l) . the infrared data have been corrected 
with a new model of extinction (J. R. D. Lepine, private 
communication). The dotted-dashed, solid and dashed 
lines show three groups of AGB stars with age of 0.25, 5 
and 9 Gyr respectively. 



Gyr and therefore the majority of them belong to the age 
range in which the AGB phase can develop significantly 
contributing to the integrated light of the cluster. Out of 
the catalog, we have selected a small sample in which the 
contribution to the total light by AGB stars is particu- 
larly strong. The photometric data of these clusters, are 
derived from the 2MASS Second Incremental Data Release 
and the companion Image Atlas, that contains about 1.9 
millions of images in J, H and K$ bands. The int egrated 
magni tudes Ji, Hi and Ki have been calculated bv lPrettd 
(2002J and kindly made available to us. Finally, in order 
to com pare our SSPs with those bv lMouhcine fc Lanconl 
( 2002]) , we consider also the integrated IR c olors for 
LMC and SMC clusters bylPerssonet all l(l98.^ used by 
iMouhcine fc Lanconl |2£)02). 

The series of Figs. El and El show the planes 
[J-H] vs [H-K], [V-K] vs [H-K], and [V - K] 
vs [J — K], respectively. In each diagram we display 
the data of IPersson et alJ l|l983j) for LMC (open cir- 
cles) and SMC ( filled c i rcles) clusters, the LMC clus- 
ters calculated bv lPrettcl 1 20021) (open squares), the SSPs 
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Fig. 15. The two color diagram [K — L] vs [J — K\. The 
sources of data and AGB models are the same as in Fig. 1141 



bvlMouhcine fc Lanconl <|2002|K the old SSPs bv lTantalol 
ill 9981). and the new ones of this study. Th e SSPs of 
iMouhcine fc Lanconl (|2002^ and lTantald l|l998t) on display 
span an age range from 100-150 Myr (the first ages at 
which the AGB contribution is significant) to 15 Gyr and 
are for the metallicities Z = 0.008, and Z = 0.02. The 
new SSPs span the same range of ages and metallicities. 

Looking at Figs . 1 1 61 1 1 71 and ITH1 we note that both clas- 
sical and new SSPs generally agree with the data, even if 
in the [J - H] vs [H - K] and [V - K) vs [H - K] dia- 
grams old SSPs span a narrower range in [H — K] and are 
too red in [ J — H] . In contrast the new SSPs have bluer 
[J — H] and [V — K] colors. Although the agreement is 
better than with the old SSPs, it is not yet satisfactory. 
An additional shift to bluer [J — H] and [H — K] is re- 
quired. This would imply that the slope of the SSP SEDs 
ought to become steeper passing from the J to the K 
pass-band th an allowed by the p r esent models. Finally 
the SSPs by IMouhcine fc Lanconl l|2002|) extend toward 
bluer colors than the other SSPs. A possible explanation 
of this unsatisfactory fit is that the theoretical models of 
M giants are not yet able to reproduce the empirical spec - 
tra of O and C stars used bv lLancon fc Mouhcind l|2002j) . 
It is worth recalling t hat the ultimate reason of t he good 
agreement achieved bv lLancon fc Mouhc ine ( 2002) models 
is the adoption of empirical spectra. More work is required 
to improve the theoretical spectra of O and C-stars to be 
included in population synthesis studies. 

Temporal evolution. Finally we look at the tem- 
poral evolution of our theoretical colors. This is shown 
in Fig. [H and H3 for the colors [J - H], [J - K] and 
[V — K]. Fig. lIHlis for the metallicity Z=0.008, whereas 
Fig-EHis for Z=0.020. The age range goes from the time 



Fig. 16. The two color diagram [J — H] vs [H — K] for 
young star clusters of the Magellanic Clouds. The open cir- 
cles ar e the LMC clusters s elected bv | Mo uhcinc fc Lanconl 
(2002) from the catalog of IPersson et al.l ljl983j) . whereas 
the filled circles are the same but for SMC clusters. The 
open squares are a few LMC clu sters whose IR colors 
have been derived bv lPrettol (2002) using the data of the 
2MASS Second Incremental Data Release and companion 
Image Atlas (see the text for more details). All the data 
have been properly reddening corrected. The lines show 
the color range spanned by SSPs of different metallicity 
and physica l input: the thin and thick dashed lines are 
the SSPs of IMouhcine fc Lanconl l|2002l) respectively for 
Z = 0.02 and Z = 0.008; the thin and thick long-dashed 
lines are our SSPs for Z = 0.02 and Z = 0.008; finally 
the thin and thick solid line s in the upper part of the di- 
agram are the old SSPs by iTantalol l|l998fr for Z = 0.02 
and Z = 0.008. The age of all SSPs is from 100-150 Myr 
to 15 Gyr. 



at which the first AGB stars are formed up to very old ages 
when the contribution of AGB stars to the integrated flux 
gets very low. For the sake of comp arison w e show also 
the colors of the cl assical SSPs by ITantalol l)l998j) and 
iGirardi et al.l l|2002j) - these latter are based on more re- 
cent stellar models calculati ons and stellar libraries - an d 
finally those of the SSPs bv lMouhcine fc Lanconl l|2002l) . 
The following remarks can be made: 

(i) The IR color s of the old SSPs bv lTantalol lll998h 
differ from those of IGirardi et alJ l|2002j) and also from 
those of the present study. 

(ii) The IR colors of the ne w SSPs marginally 
agree with the semiempirical ones bv lMouhcine fc Lanconl 
ll2002l) . In particular we note that the IR colors are in 
good mutual agreement for ages older than about lGyr, 
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Fig. 17. The two color diagram [V — K] vs [H — K] in 
the near IR. The meaning of the symbols is the same as 
in Fig. EH 
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Fig. 18. The two color diagram [V — K] vs [J — K] in 
the near IR. The meaning of the symbols is the same as 
in Fig. 



Fig. 19. The integrated colors [J — H], [J — K], and and 
[V — K] as function of the age in the range 0.1 to 20 
Gyr for SS Ps with Z=0.0 08. The dashed lines are the 
old SSPs bv lTantalol l)l998|) . The thick solid lines are the 
SSPs of this study for the same m etallicity. The dotte d 
lines are the corresponding SSPs bv lGirar di et al.l(|2002h . 
Finally, the thick long da shed line shows the SSPs of 
iMonhcine fc Lanconl j2002h for Z = 0.008. 
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whereas for younger ages, where the contribution of AGB 
stars is important, they can significantly differ. 

(iii) Finally, it is worth noticing that IR colors from 
diffe rent sources off er an embarrassing picture. Those of 
iGirardi et al.1 1 20021) are redder than any other, except for 
the [J — K) bump in lMouhcine fc Lanco rl ll2002h at about 



Fig. 20. The integrated color [J — H], [J — K], and and 
[V — K] as function of the age in the range 0.1 to 20 Gyr 
for SSPs with Z=0.02. The meaning of the symbols is the 
same of Fig. ^5] 
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Fig. 21. Upper panel: comparison between the AGB 
of various isochron es set. We plot the AGB for old 
iBertelli et alJ l|l994|) isochrones (filled, triangles). iTantakl 
lll998l) isochr o nes u sed in this work (open stars), new 
iGirardi et al.l l|2002j) isochrones (filled sq uares). Lower 
panel: comparison between [J-H] colors of [ Girardi et alJ 
20021 ) SSPs calculated wi th IGirardi et al l (1200 21) library 



of stellar spectra and with lLeimmeet" all ifl 

cal library. 



thcorcti 



l Gyr. The color s oflTantalol l)l998|) are bluer than those 
of IGirardi et alJ l|2002j) and corrected for dust effect sur- 
prisingly th ey get even bluer, w i th some excep tions. So, 
the colors of lGirardi et"afl l|2002l) . ITantald" l|l998l) and ours 
form a sequence from red to blue colors that deserves a 
careful analysis. First o f all w e ha ve to clarify the dif- 
ference between Ffantald (|l998|) and IGirardi et alJ (|2002|) 
without shells of dust. 

Spect ral libraries. A pla usible cause o f disagreement 
hetween lGirardi et alJ i|2002l) IR colors and lTantalol l|l998l) 
could be the differ e nt lib raries of stellar spectra in 
usage. IGirardi et ail l|2002|) have adopted a library in 
which some empirical spectra f or M giant star s are in- 
cluded (see IGirardi et ail 120021 for all details) . ITantald 
l|l998l) and this study adopt the ILeieune et ail (|l998l) li- 
brary in which purely theoretical spectra (even for M- 
type stars) are included. To check the effect of different 
spectral libra r ies w e have taken the isochrones/SSPs of 
IGirardi et al.1 llioO^ and recal culat ed the IR color s using 
both the IGirardi et a.lJ l)2002|) and ILeieune et all ljl99s|) 
libraries. The results are compared in the bottom panel 
of Fig. [2] limited to the case of the [J-H] color. Assigned 
the set of isochrones/SSPs . no significant differ ence arises 
passing from the libr ary of lGirardi et "all ll2002h to that of 
ILeieune et all I lj 



(1998) 

the IVassiliadis &: Wooc 
included) , but for the 



Models of AGB stars. Therefore, the cause of disagree- 
ment in the IR colors of Figs. 1191 and 1201 could be in the 
different models for the AGB phase: 

(a) Mass-loss. The first important para meter to look 
at is th e rate o f mass- loss. The SSPs of IGirardi et alJ 
120021). ITantald (Il998^) and of our study stand on 

ll993f) prescription (super-wind 
delay correction proposed by 
Vassiliadis fe Woodl l)1993|) which is not included in 
Girardi et alT l|2002). The effect of mass-loss has therefore 
little influence on the colors difference. 

(b) Path in the HRD. Another source of disagree- 
ment could be the distribution of AGB stars in the HR- 
Diagram. To this aim we compare in top panel of Fig. 1211 
the locus pred i cted by three se t s of is ochro nes/SSPs, i.e. 
IBertelli et all l|l994l) . ITantald l|l998l) and IGirardi et all 
(2002J). The comparison is limited to the age of 0.3 
Gyr and solar metallicity (Z ~ 0.02). At this age 
the contribution of the AGB to IR colors is important 
and the difference among the so urces of IR co l ors is 
stro nger. We not e that while the IGirardi et alJ J2002) 
and Bertelli et al. I lll994l) AGB is a straight line (the one 
of IGirar di et al J ll2002l) is even bl uer and f ainter than 
IBertelli et alJ |l994l)). the AGB o f ITantald lll998ft first 



coincides with IBertelli et alJ 1 19941) and eventually bends 
to cooler T e ffS. All this can be explained by the slightly 
diffe rent physics adopted in the underlying stellar mod- 
els. IGirardi et al 
loss as in lTanta" 



2002) adopt the same rate of mass- 
, ..19981) but a high e r mix ing le ngth pa- 
ramet er compared to IBertelli et al ] lll994|) and ITantald 
|l^9d). The higher mixing length given bv IGirardi et al. 
(2002) is the result of the lower, more recent opacities, 
and the calibration of their stellar m odels on the Sun (see 
IGirardi et a[|l2002l for all details). IBertelli et all l|l994l) 
and ITantald l|l99*8|) used the same stellar models and 
mixing length in turn, but adopted different expressions 
for the mass-loss rate and the relationship connecting lu- 
minosity, T e ff and c ore mass (L-TLf f-M c). The rate of 
mass-loss adopted bv IBertelli et alJ lll994l) is indeed lower 
than the one used bv lTanTald l|l998j) . This explains why 
their AGB phase extends to bright luminosities. The (L 
T e ff -Mc) relationship adopted by Bej^glh^e^jil 
stan d on the previous studi e s by IBertelli et al 
an d iGroenewegen fc de Jongl l|l993|) . whereas that used 
bv ITantald (Il998l) incorporates the TP-AGB models by 
iMarigo et al l l|l99fil) . This accounts for the different slope 
in the HR-Di a gram of t he late AGB p hase passing from 
IBertelli et all jl994l) to ITantald dl9 98l). Incidenta l ly the 
(L-T e //-Mc) relationship adopted bv IGirardi et alJ J2002) 
is much similar to that of Bertelli et al J l|l994 The differ- 
ent path in the HR-Diagram bears very much on the final 
IR colors, because cooler AGB stars should imply more 
flux at longer IR wavelengths, so that a bluer [J-H] color 
is expec t ed: in deed we have that very cool AGB stars by 
ITantald ll 19981) yield bluer [J - H] co lors, whereas the hotter 
AGB stars bv IGirardi et al.1 1|2002|) yield redder colors. 

(c) Number of AGB stars. Finally, the relative number 
of AGB in a SSP must surely play a role. To probe this 
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Fig. 22. The flux emitted by a SSP with solar metallicity 
Z = 0.02 and age of 0.3 Gyr from the stars in different 
evolutionary stages. We split the contribution to the total 
flux into five phases: main sequence (TO), from TO until 
RGB tip (T-RGB), horizontal branch (HB), asymptotic 
giant branch (AGB) and finally planetary nebulae (PN) 
that are not re presented for th e sake of clarity. The dashed 
lines are for the lTantalol lll998fl SSP, whereas the s olid lines 
are same but for the SSP of lGirardi et a l. (2002). 



effect, in Fig. [221 we plot for the SSP with solar metallic- 
ity and age of 0.3 Gyr the cumulative contribution to the 
total flux by stars in different evolutionary stages. Five 
steps are considered: up to turn-off (TO), from the turn- 
off up to the tip of the RGB (T-RGB), from this up to 
the end of the He-burning phase (HeB), and from this up 
the end of the AGB (AGB), and finally the last phase 
leading to the formation of Planetary Nebulae and White 
Dwarf (this latter is not included in our analysis owing 
to the very short lifetime). Lumping together all stages 
up to HeB we note th at the cumulative flux predicted 
bvlGirardi et al. ( 2002) is lower than that predicted by 
iTantalol l|l998|) . The situation is reversed at the final step 
when the contributi on from AGB inclusion stars i s added: 
IGirardi et al.l l|2002fl SSPs have higher fluxes than lTantalol 
( 199£$). This me ans that a higher per centage of AGB stars 
is predicted bv IGirardi et al.1 l|2002f) . The result is con- 
firmed by looking at the tabulations of the integral of the 
IMF along th e isoch rones, the so called FLUM defined by 
iBertelli et alJ l| 19941) to whom the reader should refer for 
details. The reason of it is the kind of (L-T P f f-Mc) rela- 
tionship adopted bv lGirardi et al.l l)2002j) which ultimately 
drives the evolutionary rate (relative number of stars) dur- 
ing the AGB phase. So, we can conclude that the combined 
effect of the path in the HRD an d the relative numbe r of 
AGB explain the redder colors of lGirardi et al.l 1 20021) . 



Effects of dust. Finally, there is a fundamental question 
to be clarified. Wh y the inclusion of AGB dusty shells on 
the lTantalol l|l998|) isochrones /SSPs ma kes their IR colors 
bluer than those of lGirardi et al.l l(2002h instead of making 
them redder? The explanation can be found in the effect 
of the shells of dust on the stellar radiation of AGB stars. 
Dust shifts the flux from J,H and K bands to longer wave- 
lengths, but in general its effect is stronger at the shorter 
wavelengths (J band) than at the longer ones (K band). 
More flux is driven away from the shorter wavelengths 
than the longer ones and this shift tends to increase the 
J magnitude more than the H and K magnitudes. 

The energy shift changes the flux ratio in a complicated 
way thus making it slightly higher or smaller than before. 
If the flux ratio is higher the presence of the dusty shells 
yields the effect of a redder color: for example we can 
observe as the [J-H] color tend to be bluer than one would 
expect from simple-minded considerations. For [V-K] color 
the effect is simpler because dust increases only the K 
magnitude without changing the V magnitude and thus 
producing bluer colors. 

10. Summary, discussion and conclusions 

This paper provides accurate models for the dusty shells 
surrounding AGB stars and the re-processing of the radi- 
ation emitted by the central object into the far IR. The 
aim is to derive realistic SEDs of AGB stars and SEDs of 
SSPs to be used in studies of the colors of individual late- 
type stars and the integrated colors of stellar aggregates 
going from star clusters to galaxies. 

The models for the dusty shells take into account the 
physical (density and temperature profiles and optical 
depth) and chemical structure (dust grains of different 
compositions as appropriate to the evolutionary phase un- 
der consideration) and solve the radiative transfer equa- 
tions to determine the shift of part of the coming radiation 
into the far IR region of the spectrum. Particular attention 
is payed to the transition of an AGB star from the O-rich 
to the C-rich regime taking into account recent theoretical 
models of AGB stars. 

The SEDs of AGB stars are then convolved with the 
IRAS and Johnson-Cousins broad-band pass-bands to de- 
rive the IR and far IR colors to be compared with the ob- 
servational data for individual O- and C-stars. The same 
is made for the SEDs of SSPs to be compared to the inte- 
grated colors of star clusters. Agreement between theory 
and data is good if the effect of the dusty shells is consid- 
ered. 

Finally, we compare the age dependence of the in- 
tegrated colors in the near IR for different so urces of 
SSPs based on diff e rent p hysical input, i.e. ITantalol 
l|l998l) . lo"u~ardi et all l|2002lL iMouhcine fc Lanconl l|2002F 
and ours. While the various sources more or less agree for 
the oldest ages, they differ for the youngest ages, when 
there is significative contribution from AGB stars. The 
discrepancy can be ascribed to subtle details of the mod- 
els for AGB stars proposed by different authors and to the 
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use of empirical or theoretical libraries for AGB stars. To 
conclude, more work is required to derive magnitudes and 
colors of both single stars and SSPs in the near and far 
IR. 

Extensive tabulations of spectral energies distributions 
and integrated magnitudes and colors for SSPs are avail- 
able from the authors upon request. 
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